Planar Metamaterial with Toroidal Moment 
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We observe toroidal response from a planar metamaterial comprised of asymmetric split ring res- 
onators ( ASRRs) . It is showed that a toroidal-molecule can be constructed through rational arrange- 
ment of planar ASRRs as meta-atoms via manipulating structural symmetry of the meta-atoms. 
Field maps clearly indicate that the toroidal resonance paves new electromagnetic confinement style 
in a subwavelength scale. Planar scheme of manipulating the coupling among the ASRRs may stim- 
ulate research in optical region involving toroidal multipoles. Toroidal geometry together with the 
Fano resonance made high-Q response will have enormous potential applications in low-threshold 
lasing, cavity quantum electrodynamics and nonlinear processing. 
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Toroidal moment is of fundamental interest in the 
physics community [1 . The toroidal response produced 
by currents flowing on the surface of a torus along its 
meridians is firstly considered by Zel'dovich to explain 
parity violation in the weak interaction [2] . The so called 
anapoles [2] have gained acceptance in nuclear and par- 
ticle physics [31 3] . Theoretical calculation has also pre- 
dicted toroidal dipole in a kind of molecular structure 
[5]. Toroidal multipoles are different from electric or 
magnetic multipoles in traditional multipole expansions. 
Particularly, non-radiating feature of the toroidal geom- 
etry is quite unique which comes from its fantastical field 
localization configuration in a head-to-tail manner. How- 
ever, the scattering of external fields from toroidal multi- 
poles in nature is often much weaker for observation. So, 
suppressing traditional electric and magnetic multipoles 
and enhancing toroidal response by rationally manipulat- 
ing structural symmetry are basic and essential towards 
a toroidal geometry. 

Since there exist profound analog between the weak 
and the electromagnetic interactions, if we can recon- 
struct this unusual interaction in manipulating electro- 
magnetic wave which is important in our modern life? 
Metamaterial is just a concept of artificially construct- 
ing micro-structured materials with fancied properties 
unattainable in nature, these artificial materials usually 
consist of subwavclcngth-sized metallic resonant building 
blocks as meta-atoms [Si [7], which yield electric or/and 
magnetic responses. The local resonances of the meta- 
atoms and mutual coupling among them enable the ma- 
nipulation of electromagnetic waves at a subwavelength 
scale. 

Toroidal metamaterial was first theoretically proposed 
by Marinov et al. in 2007 in investigating a 3D-array 
of toroidal solenoids [5], later Papasimakis et al. ex- 
perimentally studied a toroidal wiring consisting of four 
rectangular loops and analyzed its response in terms of 
multipole moments [9 , however toroidal moment is not 



dominant in the structure. It was found toroidal meta- 
materials can make negative index of refraction and ro- 
tate polarization state of light 8,9 . In 2010 toroidal 
dominated response in microwave region was first experi- 
mentally demonstrated by three-dimensionally arranging 
four split ring resonators (SRRs) in a unit cell of toroidal 
symmetry [10 . Although there are some attempts in 
manufacturing three-dimensional SRR structures, it is 
of challenge not only at optical frequencies but also at 
microwave frequencies [10] . Stacking of two-dimensional 
planar plasmonic metamaterial structures is assumed to 
be the pathway towards the third dimension in the op- 
tical regime. These metamaterials have exhibited exotic 
properties and potential applications: artificial magnetic 
composite [IT] , perfect absorber [12] , electromagnetically 
induced transparency medium [13-16 . Planar structure 
based scheme is not limited at microwave band, but also 
show performance at Tera-hertz band [17j and even op- 
tical range [TgH20] . 

In this paper, we present a planar scheme for toroidal 
metamaterial, with the unit cell comprised of four asym- 
metric split ring resonators (ASRRs), we show that a 




FIG. 1. (a) Schematic of the unit cell for the ASRR-based pla- 
nar toroidal metamaterial and corresponding electromagnetic 
excitation configuration (with polarization direction along x 
axis), (b) Photograph of the toroidal metamaterial slab, geo- 
metric parameters are denoted by black letters in Fig. 1(a). 
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toroidal-molecule can be constructed through rational ar- 
rangement of planar ASRRs as meta-atoms via manipu- 
lating structural symmetry among the meta-atoms. Field 
maps clearly indicate that the toroidal resonance paves 
fantastical electromagnetic field confinement style in a 
subwavelength-sized head-to-tail channel of the toroidal 
geometry. Planar scheme by manipulating the coupling 
among the ASRRs may stimulate research in optical 
region involving toroidal multipoles. Toroidal geome- 
try together with the Fano resonance made high-Q re- 
sponse will have enormous potential applications in low- 
threshold lasing, cavity quantum electrodynamics and 
nonlinear processing. 

The designed planar structure is comprised of two 
metallic layers and a dielectric spacer layer. A unit 
cell of our model slab or the toroidal-molecule is shown 
stereographically in Fig. [ija), the toroidal- molecule is 
constructed through rational arrangement of four asym- 
metric split square ring resonators (ASRRs) as meta- 
atoms via manipulating structural symmetry of the meta- 
atoms. Pairs of ASRRs in the same layers are of mirror 
symmetry about the xoz plane (er) ; pairs of ASRRs in the 
different layers are of 180° rotational symmetry about the 
y axis (C2). We note that throughout our study meta- 
materials are illuminated by x-polarized electromagnetic 
waves both in experiments and calculations as illustrated 
in Fig. pa). 

Toroidal-molecule in our model system is constructed 
from four contiguous ASRRs as meta-atoms, since these 
meta-atoms are in close proximity, the nearest-neighbor 
mutual interactions or hybridizations [2TH24] between 
them play determinative role in their collaborative re- 
sponse to external excitations. We firstly investigate re- 
sponse of single ASRR and then mutual coupling inter- 
actions in pairs of the ASRRs. We performed all numer- 
ical calculations with a brute-force finite-difference-in- 
time-domain (FDTD) [25] electromagnetic solver (Gal- 
lop) [55]. For a metamaterial with its unit cell only one 
ASRR, it is showed that the trapped mode can be excited 
even under normal electromagnetic incidence condition 
in a recent study [57] , the trapped mode originates from 
symmetry breaking to some extent of the metallic ele- 
ment of a metamaterial. Small structural asymmetry in- 
troducing makes a degeneration from dark mode |28j to 
the trapped mode, or we can say weakly coupled nature 
to external excitation of ASRR is an inherited attribute 
from the dark mode. For the single ASRR made unit cell 
case, frequency plot performance of the transmission is 
shown as black solid line in Fig. [2j The dip position near 
10 GHz indicate the trapped mode of ASRR, the resonant 
response appear to have Fano line-shape [2"5H5T] . result 
in destructive interference between trapped sharp reso- 
nance and free space propagating continuum-like spec- 
trum. Meanwhile, we noted that the single ASRR is 
of weak excitation, again verify that the trapped mode 
of ASRR is weakly coupled to external excitation since 



electromagnetic wave can easily leak through the unpat- 
terned area of the meta-surface. 

Then we consider coupling interactions between the 
ASRRs. Numerical simulations were first carried out for 
the structure of in-plane ASRR pair [i.e., 5ft and a(5ft) 
in Fig. [lja)] as fundamental element, the transmittance 
spectrum is represented by the blue solid curve in Fig. [2] 
in which only one resonance is clearly observable. We 
call this kind of interaction as horizontal coupling, the 
single peak response comes from in-plane mirror symme- 
try, only mirror symmetric surface current can be excited 
under normal incidence because anti-symmetry surface 
current can not be excited since there is no phase re- 
tardation, the anti-symmetric mode is now a dark mode 

For the structure of stacked ASRR pair [i.e., 5ft and 
6*2(3?) in Fig. [lja)] as fundamental element, the trans- 
mittance spectrum is represented by the red solid curve 
in Fig. [2] there now exist two resonances. We call this 
kind of interaction as vertical coupling, the two reso- 
nances are far away from the original magnetic trapped 
mode of ASRR, which is due to extreme overlap of the 
two magnetic trapped modes of ASRRs in the vertical 
coupling configuration. We also note that for the low- 
frequency mode magnetic oscillations of the ASRRs are 
in-phase, this is crucial for toroidal response in our pla- 
nar metamaterial system, the 2-fold rotational symmetry 
of the mcta-molecule made the in-phase response at low 
frequency. 

We construct a meta-molecule involving previously 
mentioned horizontal coupling and vertical coupling with 
four ASRRs for toroidal response. Figure [3] presents the 
FDTD calculated and experimentally measured trans- 
mission spectra. The calculated intensity of transmit- 
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FIG. 2. Comparison of FDTD calculated transmission spec- 
tra for planar metamaterial unit cell with only single ASRR 
(black curve), vertical coupled ASRR-pair (red curve) and 
horizontal coupled ASRR-pair (blue curve). 
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FIG. 3. Experimentally measured (black curve) and FDTD 
calculated (red curve) dispersion of transmitted power of the 
designed planar toroidal metamaterial. 

ted wave (red solid line) shows clearly two Fano shaped 
resonances come from destructive interference between 
collective sharp response of the meta-molecule and free 
space propagating continuum-like spectrum. 

The two resonances have their characteristic frequency 
nearly the same as vertical coupled ASRR pair (red solid 
curve in Fig. j2J), frequency of the upper resonance at 
12.905 GHz (marked with the peak frequency) is deter- 
mined by the out-phase vertical coupling, which is clearly 
observable as in the color map of z-component magnetic 
field in Fig.[4]Jb), the coupling in this manner causes a net 
nonzero magnetic component along y axis; meanwhile, 
frequency of the lower resonance at 6.378 GHz (marked 
with the peak frequency) is determined by the in-phase 
vertical coupling, in-phase coupled magnetic mode [see 
color map of the z-component magnetic field in Fig.|4^a)] 
together with horizontal coupling under normal electro- 
magnetic incidence made a circumfluent magnetic field 
or magnetic votex, this is rightly a picture of the toroidal 
dipolar [10] which shows a head-to-tail magnetic config- 
uration. 

We fabricated a sample with exactly the same param- 
eters as our theoretical model. A printed circuit board 
(PCB) with metal sheets (0.035 mm-thick) attached to 
both sides of a 0.787 mm-thick commercial Taconic TLX- 
8 dielectric substrate with the permittivity e r = 2.55 is 
employed. The meta-molecule array has a lattice con- 
stant of P x = 5 mm and P y — 10 mm. The metal strip 
width is w = 0.5 mm wide, the outer width of the ASRRs 
are L = 4 mm, inner gap between the two in-plane AS- 
RRs is g — 0.8 mm, and a s = 1 mm wide split is opened 
on the original square ring with a distance Z = 1.9 mm 
away from the center sideline, overall lateral size of our 
planar metamaterial slab is 540 mm x 440 mm. The 
transmission spectra through our sample were measured 




FIG. 4. Calculated spatial distributions of vectorial magnetic 
field for correspondingly toroidal mode and magnetic mode 
plotted over the color maps of z-component magnetic field 
at peak frequencies 6.378GHz (a) and 12.905GHz (b) respec- 
tively. 



with a vector network analyzer Agilent 8722ES in an ane- 
choic chamber. 

The measured results (black solid curve in Fig. [3]) meet 
quite well with the theoretical predictions. Two reso- 
nances of our model system are all clearly observable in 
the experimental spectrum. However we found a curi- 
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FIG. 5. Frequency dependence of radiating power for multi- 
poles from the planar toroidal metamaterial: toroidal dipole 
(red curve), magnetic dipole (blue curve) and electric dipole 
(cyan curve). The inset illustrates radiating strength of the 
toroidal dipole and electric dipole around 6.378GHz. 



ous and abnormal phenomenon: the low-Q coupled mag- 
netic resonance at high frequency is less obvious while 
the high-Q toroidal resonance at low frequency is quite 
obvious. We believe that this is due to the exotic field 
localization configuration of the toroidal geometry. This 
can be checked by looking at the evolution of magnetic 
field inside the structure as calculating spatial field dis- 
tribution showed in Fig.|4j in the toroidal mode (at 6.378 
GHz) case, most field is restricted to flowing inside the 
head-to-tail channel while for the coupled magnetic mode 
little field is confined inside the channel, corresponding 
to in-phase and out-phase vertical coupling as illustrated 
in the color maps of z-component magnetic field which 
show continuum or reversion. It can be see that strongly 
enhanced field is adhered and localized on metal surfaces 
in the coupled magnetic mode (at 12.905 GHz) as in pre- 
vious metamaterial structures, Planar toroidal metama- 
terial offers strikingly different style for electromagnetic 
is of convincing manifestation about the abnormal low 
loss phenomenon of the toroidal mode. The findings have 
beneficial for low Joule loss metallic photonic mctamate- 
rials |52] . 

Radiating power for induced multipole moments is cal- 
culated through integrating spatial distributed fields over 
volume of the planar meta-molecule with the current 
density formalization [TU] . Figure [5] shows the radiat- 
ing power of three dominativc multipolcs of the planar 
metamaterial system: electric dipole (cyan curve), mag- 
netic dipole (blue curve) and toroidal dipole (red curve) . 
As can be seen from Fig. [5] electric dipole P x shows large 
components in the full frequency band, which reveal the 
fact that our planar metamaterial system is electric ex- 
cited and the planar structure is hard to coupling with 



the incoming parallel magnetic component. Generally we 
notice that magnetic dipole M y and toroidal dipole T x 
show considerable scattering contributions at 6.356 GHz 
and 12.930 GHz which are just near the two transmission 
peak frequencies. Though ^-component of electric dipole 
moment P x show strongest contribution at the 12.905 
GHz transmission peak frequency, the y-component of 
the magnetic dipole moment M y also provide significant 
contribution which can be checked in the time-domain 
simulated field patterns. While for the 6.378 GHz trans- 
mission peak, the x-component of the toroidal dipole mo- 
ment T x plays decisive role in comparison with electric 
dipole and magnetic dipole, inset in Fig. [5] clearly illus- 
trates radiating strength of the toroidal dipole moment 
T x and electric dipole moment P x around 6.378 GHz, we 
see that at 6.378 GHz radiating power for the toroidal 
dipole moment is three orders more of magnitude larger 
than that for electric dipole moment, this is due to elec- 
tric dipole radiation shows a Fano line-shaped spectrum 
that comes from radiating interference of induced current 
in the trapped mode. All these indicate that the exper- 
imentally observed transmission peak around 6.378 GHz 
is the toroidal moment response of our planar metama- 
terial. 



In summary a planar metamaterial structure com- 
prised of ASRRs is proposed for observation of toroidal 
moment response. It is showed that a toroidal-molecule 
can be constructed through rational arrangement of pla- 
nar ASRRs as meta-atoms via manipulating structural 
symmetry of the meta-atoms. The characteristic spec- 
trum is experimentally verified at microwave range. Field 
maps and multipole moment calculations clearly indicate 
that the toroidal resonance paves new electromagnetic 
confinement style in a subwavelength scale. The pla- 
nar toroidal structure can be generalized to terahertz, 
infrared, and visible frequencies particularly for low Joule 
loss metallic photonic metamaterials, the novel field lo- 
calization configuration with a high-Q response and low 
loss have beneficial for myriad photonic applications, 
e.g., low-threshold lasing, cavity quantum electrodynam- 
ics and nonlinear processing. 
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